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Abstract — In cellular wireless networks, user association refers 
to the problem of assigning mobile users to base station cells 
- a critical, but challenging, problem in many emerging small 
cell and heterogeneous networks. This paper considers a gen- 
eral class of utility maximization problems for joint optimiza- 
tion of mobile user associations and bandwidth and power 
allocations. The formulation can incorporate a large class of 
network topologies, interference models, SNR-to-rate mappings 
and network constraints. In addition, the model can applied 
in carrier aggregation scenarios where mobiles can be served 
by multiple cells simultaneously. While the problem is non- 
convex, our main contribution shows that the optimization admits 
a separable dual decomposition. This property enables fast 
computation of upper bounds on the utility as well as an efficient, 
distributed implementation for approximate local optimization 
via augmented Lagrangian techniques. Simulations are presented 
in heterogeneous networks with mixtures of macro and picocells. 
We demonstrate significant value of the proposed methods in 
scenarios with variable backhaul capacity in the femtocell links 
and in cases where the user density is sufficiently low that lightly- 
used cells can reduce power. 

Index Terms — cellular wireless systems, user association, dual 
decomposition, inter-cell interference coordination, carrier ag- 
gregation, 3GPP LTE. 



I. Introduction 

In cellular wireless networks, user association refers to the 
problem of assigning each mobile terminal a serving base 
station cell. In traditional cellular systems, user association is 
relatively straightforward in that mobiles can simply connect 
to the cell with the strongest received signal strength. This 
simple selection policy maximizes the SNR to each user and 
hence the rate per unit bandwidth. 

However, in many emerging cellular network deployment 
models, user association requires consideration of other factors 
in addition to signal strength. For example, to scale capacity 
in a cost effective manner, traditional macrocellular networks 
are being increasingly supplemented with overlays of smaller 
micro- and picocells (TJ, (2). The resulting heterogeneous 
networks (HetNets) will consist of cells with vastly differ- 
ent sizes and require load balancing to encourage mobiles 
to connect to smaller cells that would otherwise be under- 
utilized (3), (4). In addition, these HetNets may also contain 
open femtocells pi, (6) where the wired connectivity would 
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be provided by third parties other than the cellular provider. 
Third party backhaul capacity may be highly variable and may 
also need to be considered in the user association decisions. 
Finally, the state-of-the-art cellular standards offer a number of 
advanced intercellular interference coordination (ICIC) mech- 
anisms including subband scheduling and beamforming (7). 
As a result, the throughput a user can experience in any one 
cell can depend in a complex manner on the particular choice 
of user assignments and resource allocation decisions in other 
cells. 

To understand these issues, this paper considers a general 
class of problems for joint optimization of user association and 
bandwidth and power allocations for interference coordination 
in cellular wireless networks. The joint user association- 
interference problem is formulated in the classical framework 
of utility maximization (8), (9). Specifically, the throughput 
to the users is regarded as a function of the cells that they 
are assigned to as well as the resources (power or bandwidth) 
allocated to the users in the cells. The bandwidth and power 
allocations in one cell may affect the interference levels in 
other cells. The problem is to match the users to cells and 
allocate resources within each cell to control the inter-cellular 
interference and maximize some system- wide utility function 
of the rates. 

Two variants of the problem are considered: (a) Multiflow 
optimization where mobiles can be served by multiple cells 
simultaneously as in the case of 3 GPP LTE with carrier aggre- 
gation (CA) (10), (TTJ; and (b) user association optimization 
where each mobile can only be served by one base station 
cell at a time as in the case of standard 3 GPP LTE without 
advanced CA capability. 

Both problems are, in general, non-convex. However, while 
we cannot find an optimal solution, our main contribution 
shows that under very general formulations, the multiflow 
optimization admits a separable dual decomposition. This 
separability property has several key implications: First, for 
any given values of the Lagrange parameters, the dual problem 
can be solved easily and hence one can minimize a dual upper 
bound exactly. Therefore, we can obtain a computable upper 
bound against which we can compare the performance of 
any practical, suboptimal algorithm. Secondly, although the 
problem may have a non-zero duality gap (i.e. the problem 
may not be strongly dual), the separability of the dual prob- 
lem enables efficient implementation of various augmented 
Lagrangian techniques for constrained maximization. Finally, 
dual techniques have the benefits that they readily lead to 
distributed implementations where the base stations broadcast 



the dual parameters and the mobiles update the cell selections. 
In the multiflow problem, the Lagrange parameters have a 
natural interpretation as access prices. 

For user association, we recognize that the problem is 
essentially classic route selection in networking theory. Route 
selection is known in general to be NP-hard, but we follow a 
standard procedure where we use the solution to the multiflow 
problem as an upper bound and then truncate the solution for 
the single path constraint 1 12]— 1 14]. In this way, we can apply 
dual decomposition plus one additional truncation step for the 
user association problem. This mutliflow approximation has 
also been used in cellular user association in (4) (see below). 

An appealing feature of our approach is that the dual 
separability applies under very general circumstances. Specif- 
ically, our formulation incorporates a general linear mixing 
interference model, first presented in [15] for ICIC problems 
in femtocells. Although we will only study bandwidth and 
power allocations in this paper, the model can incorporate 
- with some extensions - other ICIC techniques including 
subband allocations and beamforming. In addition, our only 
restriction on the utility functions is that they are separable. 
The dual decomposition will hold even under non-concave 
utilities. Finally, unlike previous formulations, the approach 
here can incorporate network constraints (provided they are 
known) including limits in the backhaul - an important feature 
for third-party offload. 

A. Related Work 

Early work in the area of optimal user association includes 
the classic paper by Hanly [16] that found a simple iterative 
algorithm for jointly optimal user association and uplink power 
control in CDMA systems. Downlink versions of this problem 
were later considered in (17). Although we will focus on 
bandwidth allocation for OFDMA systems in this paper, our 
methods apply to general linear mixing interference models, 
and thus power control may also be in considered in the future. 

The recent interest in user association is however not for 
power control, but rather optimized load balancing in het- 
erogeneous networks as described in (3), (4), (18). Solutions 
in 3GPP standards bodies have focussed on simple "range 
expansion" techniques that apply a fixed bias in the cell 
selection procedure to shift mobiles to the smaller cells (19), 
(20). These techniques can be supplemented with adaptive bias 
based on SINR [3] and iterative transmit power selection (18). 
The cell-site selection problem is addressed as utility-based 
optimization problems with backhaul capacity constraints 
in |2T) , while (22) proposes a distributed algorithm based 
on an optimal load of base stations. Although our evaluation 
methodology is based on standard simulation models such 
as (23) , stochastic analyses have been provided for these range 
expansion techniques in (24). The methodology in this paper 
follows most closely to the work [4] that used a similar 
optimization technique with a multiflow upper bound and dual 
decomposition. The work |4] however does not consider joint 
optimization of user association with power and bandwidth 
control. 

One of the interesting findings in (3), (4), (T8) was that static 
range expansion works remarkably well and more sophisti- 



cated optimization approaches offer little benefit. However, 
here we will see that, when combined with joint bandwidth / 
power optimization, or in the presence of network constraints, 
optimized user association can offer larger gains in some 
scenarios. For example, as our simulations will indicate, when 
the network loading per cell is light or when the backhaul 
capacity is limited, unused cells can reduce power improving 
the interference conditions in other cells. 

To perform the joint optimization of user association and 
interference coordination, our work synthesizes several well- 
known network optimization methods. The use of multipath 
upper bounds for optimized route selection is a classic tech- 
nique in networking theory that dates at least back to (12), 
[13 1 - see the text [14] for a historical survey as well as some 
testable conditions under which the multipath upper bound 
is provable optimal. See also (8), (25)-(27). Joint power and 
congestion problems is given in the well-known paper [28] . 
See also (29) and the monograph (30). This works also uses 
dual decomposition techniques, but applies to a potentially 
more general linear mixing interference model that includes 
power control as a special. However, while more general, our 
methodology has no guarantee of global optimality. 

Utility optimization approaches to the multipath problems 
have been extensively studied including both primal and 
dual solutions with close connections to rate control (8) and 
network pricing [25]. The dual decomposition methods used 
in this paper are closest to (26), (27) which explore augmented 
Lagrangian and method of multiplier techniques. 

II. General Flow and Interference Control 
Problem 

Before describing the user association and multiflow prob- 
lems in cellular networks, we first describe a general joint 
flow control and interference optimization problem. We will 
apply this general formulation to cellular networks in the next 
section. 

For the general problem, we consider a system with S flows 
with rates r 39 s = 1, . . . , S where a subset of the flows are 
constrained by interfering wireless links. Index the wireless 
links by £ = 1, . . . , L and let T(f) denote the set of flows 
that share the wireless link L We assume that wireless link 
capacity places a constraint on the rates through the link of 
the form 



ser(£) 



r s < C£(x£,ze), £ 



(D 



where X£ denotes a resource allocation variable for the link, 
Z£ is the interference level on that link and Cg(xg^zg) is the 
wireless link capacity. Following [15], our key assumption will 
be that the interference level on the links are a linear function 
of the resource allocation variables on the other links so that 
we can write 

z = Gx, (2) 

for some gain matrix G. In addition, we assume that there are 
some linear constraints on the resource allocation variables of 
the form 

A x x < b x , (3) 
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Wireless link b 
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Fig. 1: Simple example of a network with three flows and 
two interfering wireless links for joint power control and rate 
optimization. 



for some matrix A x and vector b^. In the sequel, it will be 
convenient to rewrite the L constraints ([]]) in a single matrix 
form 

Er<C(x,z), (4) 

where E is the incidence matrix with components 



_ J 1 if seT(£), 



Eis ~ "' ifsf£T(£), 
and C(x, z) is the capacity vector 

C(x,z) = (CiOi,zi),--- ,C L (x L ,z L )) T . 

In the user selection problem below, xg will be a quantity 
proportional to the power allocated to a wireless link L Hence 
the interference powers zg are naturally described as a linear 
function of the form ([2]) where the matrix G would depend 
on the path losses between the transmitters and receivers. The 
link capacity ([T]) could then be a function of the signal-to- 
noise ratio xg/zg. As described in [15], the linear mixing 
interference model ^ can also be used for studying more 
complex interference scenarios including links with subbands 
or beamforming. 

As discussed in the Introduction, we also wish to incor- 
porate rate constraints in the wired network. To this end, we 
assume the network constraints are of the form p4| 



A r r < b r 



(5) 



for some matrix A r and vector b r . 

Our goal is to select a resource allocation vector x and flow 
rate vector r to maximize some utility function of the form 



U(r) 



S 



U a (r a ), 



(6) 



where U s (r s ) is the utility of the rate on flow s. Our 
formulation will not require that the utility functions are 
concave, only that they are continuous and monotonically 
increasing in r. 

To illustrate the formulation, Fig. [T] shows a simple example 
for joint power control and rate optimization. The network has 
S = 3 flows with rates r\, r 2 and r$. Two of the links, labeled 
links a and b, are wireless links that interfere with one another. 
The transmit signal powers on the links are denoted x a and 



x b , while z a and z b represent the interference powers. In this 
case, 

where G a b and G ba are the cross gains between the transmit- 
ters and receivers. The SINRs on the two links would be given 

by 

G aa x a G bb x b 

la = ; , 76 = : , 

z a + n a z b + n b 

where n a and n b are the noise powers and G aa and G bb are 
the gains on the intended links. If the link capacities C a and 
C b can be modeled as functions of the SINRs j a and 75, then 
we can write 

r\ < C a = C a (x a , z a ), r 2 + r 3 <C b = C b (x b , z h ), 

which fits in the format of ([I]). Under this model, one can then 
attempt to jointly attempt to optimize the transmit powers x a 
and x b along with the flow rates 7*1, r 2 and r% to maximize 
some utility. 

To summarize the optimization, let 6 denote the vector of 
all the decisions variables 



= (r,x,z). 
and let 6 denote the maximum values 
= (r,x,z). 



(7) 



(8) 



With some abuse of notation, we let U(0) = U(r), so that we 
can rewrite the optimization as 



max_£/(0) s.t. g(0) < 0, 

0<6<6 



(9) 



where g{0) is the vector of constraints 

9(0) ■ 
and is some upper bound on the variables. 



£*-rpL KJrp 

Gx-z 

Er-C(x,z) 



(10) 



III. User Association and Multiflow in Cellular 
Networks 

Having described the general flow and interference control 
problem, we now turn to special case of user association 
and bandwidth optimization in cellular systems. Our model 
considers the downlink of a cellular network with Nms mobile 
stations, MS z, i = 1, . . . , A^ms and AT BS base stations, BS j, 
j — 1, . . . , A^bs- We suppose that there are a total of L flows, 
with each flow £ passing through exactly one BS-MS pair. 
Although a single base station may support multiple flows, 
we will logically regard each flow £ as a separate wireless 
link with rate 77 and bandwidth allocation xg. For simplicity, 
we will assume only the BS-MS links are wireless with 
interference constraints. Any other (wired) links are assumed 
to have fixed capacity constraints. We let Tms(^) denote the 
set of flows received by MS i and TbsO) t> e me set °f fl° ws 
transmitted by BS j. 

To illustrate the model and notation, Fig. [2] shows a simple 
example of a network with two mobile stations, MSI and MS2, 



Candidate base 
stations (macros, 
picos, femtos) 




Fig. 2: Simple example network where two mobiles, MSI and 
MS2, must each select between two base stations. Since each 
of the two MS's have two candidate flows, there are a total of 
L = 4 flows with sets T M s(l) = {1,2}, T M s(2) = {3,4}, 
Tbs{1) = {1}> Tbs(2) = {2,3} and F B s(3) = {4}. 



in a typical cellular configuration. Mobiles MSI and MS2 
received data from sources SI and S2 in the public Internet. 
Following a standard cellular architecture such as 3 GPP LTE 
(3TJ, the data is assumed to arrive from the public Internet to 
a common gateway which then tunnels the traffic through an 
operator-controlled core network to the base station selected to 
the serve the mobiles. In the simple example in Fig. [2] we have 
three candidate base stations shared amongst the two mobiles. 
Since each of the two MS's have two candidate BSs, there are 
a total of L = 4 flows, each flow following a path of the form: 
source — » gateway — >• BS — )> MS. 

Now returning to the general case, let fi be the total rate to 
MSz, which must satisfy the constraint 



n < 



E 

eer MS (i) 



re, 



1,...,JVms- 



(11) 



We let r be the vector containing both the total and individual 
flow rates 



(ri 



^N MS ,r u ...,r L ) 



(12) 



We assume the utility is a separable function of the total rates 
to each MS so that 



U(t) 



Nms 

E 

i=l 



Ui(fi), 



(13) 



for some utility functions Ui(fi). 

To model the bandwidth constraints and interference, we 
assume that each BS j has a total available bandwidth Wj and 
radiates a fixed power per unit bandwidth Pj. We let xg be 
the bandwidth allocated to the flow on link £. The bandwidth 
allocations at BS j must thus satisfy the constraint 

J2 x £ <Wj, j = l,...,7V BS . (14) 

ter BS (j) 

Also, the total power radiated by BS j is given by 

Pn 



3 Yl X£ ' 



Now suppose that a flow £ is served by BS j and received by 
MS i. Then, the total interference power received on that flow 
is given by 

z £ = ^2H £k P k ^2 x ^ ( 15 ) 

Mi *er BS (fc) 

where H^ is the channel gain from BS/c to the receiver for 
flow £. Since flow £ is served by BSj, the SINR on the link 
is then 



li{zt) 



zt + Nt 



where Ne is the thermal noise at the receiver. We can then 
assume the capacity constraint on the £-\h link is of the form 



r £ < C(x£,zi) := xipi(^i(zi)), 



(16) 



where pg is the spectral efficiency (rate per unit bandwidth) on 
the £-\h link as a function of the SINR. The spectral efficiency 
pi is multiplied by the allocated bandwidth xg to yield the link 
capacity Cgix^zg) in ( [T6| ). Our formulation will permit any 
continuous spectral efficiency function. 

Under this model, we consider two problems: 

• Multiflow optimization: In this case, there are no ad- 
ditional restrictions on the rate vector r so that the 
mobile can be served from multiple base station cells 
simultaneously. In 3 GPP LTE- Advanced, this scenario 
would correspond to carrier aggregation |T0) , (TTJ. 

• User association optimization: In this case, we as- 
sume each mobile can be served by only one cell at a 
time, which is the standard model for LTE without the 
advanced carrier aggregation feature. This requirement 
can be modeled via an additional constraint: for all 
i = 1,...,JVms, 



ri = for all but one i G TmsW- 



(17) 



^r BS (j) 



Under this constraint, the selection of the index £ such 
that r£ 7^ 0, determines which BS will serve MS i. 

Now, the multiflow optimization is exactly a special case of 
the general flow and interference control problem in Section [H] 
Specifically, let x G IR L be the vector of bandwidth allocations 
and z G R L be the vector of interference powers. The band- 
width constraints ( fT4| ) can be written of the form A x x < b x 
for some A x and b^. The network constraints will be written 
in the form A r r < b r , which can include constraints on either 
the flow rates or total rates. The interference powers fl5] ) can 
be written as z = Gx for some gain matrix G and the capacity 
fT6] ) can be written in the form ([4}. 

On the other hand, the single path constraint ^Tf\ for 
the user association optimization problem cannot be directly 
placed into the form of flow and interference control problem 
in Section [HI As we will discuss in the next section, we will 
use the multiflow optimization as an upper bound for this 
problem and then truncate the output to obtain a feasible single 
path solution. 



IV. Approximate Optimization via Dual 
Decomposition 

A. General Flow and Interference Control 

We now describe potential methods for solving the opti- 
mization problems. We begin with the general flow and in- 
terference control problem in Section [TT| The maximization of 
([9]) in this problem is non-convex in general since the capacity 
function C(x, z) may have an arbitrary form. In addition, the 
utility function U(y) may be non-concave. However, since 
both the objective functions and constraints in ^ admit a 
separable structure, it is natural to attempt to solve the ([9]) via 
dual decomposition. 

To this end, corresponding to the constrained optimization 
([9]), define the Lagrangian [32] 



L{0,n):=U{e)-ti T g{0) 



= U(t) + n xl (h x - A x x) + n rl (b r - A r r) 
+ n zT (z - Gx) + n cT (C(x, z) - Er) , 



(18) 



where \i > are the dual parameters which are partitioned 
conformably with g(0) in < |T0] > as 



n = (n x ,n r ,n z ,n c ) 1 



(19) 



The following simple lemma provides our main motivation for 
considering duality -based optimizations: namely that the dual 
maximization is separable: 

Lemma 1: Let $(0) be any separable function of the form 

*(0) = $(r,x,z) := X><fo) +$>! (xi,zi), (20) 

i i 

for some functions </>[(•) and </>f (•)• Then, the maximization 
of the Lagrangian L(0, /x) in ([18]) augmented by $(0), 

0(/x) = (r(/x),x(/x),z(/x)) := argmax L(0,/x) - $(0) 

0<0 L J 

(21) 

has components given by the solutions to the optimizations 



r s G argmax 



C/ S (r s ) - #(r a ) - (E T /x c + A^ r ) s r £ 



s.t. < r s < r s 



(22) 



and 



(£*,2i) G argmax il\Ci(xi,zi) - 4%(xi) - 4>\{zi) 

-(A^/x* + G T /ji z )iXi + /if ^ 
s.t. < x^ < x^, < ^ < zg. 



(23) 



Proof: This result follows immediately from the separable 
structure of the objective function and constraints. ■ 

Thus, the vector- valued optimization ( |2T] ) decomposes into 
5 one-dimensional optimizations over the variables r s and L 
two-dimensional optimizations over the variables {x^zij. In 
many cases, we will see that these have explicit closed-form 
solutions or can be easily computed via a simple line search. 

There are two key consequences to this separability result: 



a) Computation of an Upper Bound: First, the well- 
known weak duality property (32) states that an upper bound 
to the maximum of the utility is given by 



max U(0) < minL(u), 

6 : g(0)<0 /x>0 



where L(fi) is the dual maxima, 



L(fjb) := maxL(0,/Li). 
e<e 



(24) 



(25) 



= arg max 

9<G 


' L{e, ii*)- 


-$(0,0*) 


(26a) 


= ^-ag{6 t+l ), 




(26b) 



This bound applies even when U{6) is non-convex. Moreover, 
L(/jl) is also always convex in fi with subgradient g(0(ii)). 
Since Lemma [l] shows that the maxima 0(/x) and L(fi) can 
be easily evaluated for each fj,, the upper bound in ( [24] ) can be 
minimized via convex programming, thereby providing a com- 
putationally tractable upper bound on the utility. Of course, 
since U(0) and/or the constraints may be non-convex, the 
bound may not be tight (i.e. the problem is not strongly dual). 
Also, in absence of the augmenting term, the subgradient may 
be non- smooth, meaning the convex optimization may have a 
slower (although still feasible) convergence rate (33), (34). 

b) Augmented Lagrangian Optimization: Secondly, since 
the maximization ( |2T] ) can be computed for any fi and separa- 
ble augmenting function <£, one can employ a wide range of 
augmented Lagrangian optimization methods for finding local 
optima for the constrained maxima problem ([9]). Augmented 
Lagrangian methods |32| generate a sequence of estimates f 
typically through the iterations of the form 



where a > is a step-size and ^(O^O 1 ) > is a so-called 
augmenting function that smooths the optimization. Separable 
forms of the augmenting function can include the quadratic 

$(0,0*) =e||0- 0*||| for some e > 0. 

B. Multiflow Optimization 

Since the multiflow cellular optimization in Section [TIT] is a 
special case of the general flow and interference control prob- 
lem, we can directly apply the results above in Section |IV-A| 
to compute upper bounds on the utility as well as efficiently 
implement augmented Lagrangian techniques. 

In addition, in the cellular context augmented Lagrangian 
iterations of the form ( [26] ) also admit a simple distributed 
implementation - a feature also found in the similar user 
association optimization algorithm in (4), as well as many 
duality-based power control algorithms (see the monograph 
(30)). Specifically, it can be verified that, using the separable 
structure of the optimization in Lemma [T] the rates th and 
bandwidth allocations xg and target interference level zg for 
each link can be computed at the base station serving that 
flow. The base station only needs local information on the 
gain matrix G and dual parameters. The base stations can 
then update the dual parameters fi as in ( |26b| ) and broadcast 
the values to the neighboring base stations potentially through 
the mobiles. 



TABLE I: Simulation parameters 



Max-RSSI vs. Range Expansion vs. Optimizer (P=0) 



Parameter 


Value 


Macro layout 


19 hexagonal cell sites with wraparound, 3 cells 
per site. Inter-site distance=500 m. 


Pico and UE distri- 
bution 


10 picos and 10 or 25 UEs per macrocell, uni- 
formly distributed. 


TX power 


46 dBm (macro), 30 dBm (pico) 


Path loss 


macro: 131.1 + 42.8 log 10 (#), 
pico: 145.4 + 37.5 log 10 (ifc) 


Carrier frequency 


2.1 GHz 


Shadowing 


8 dB lognormal std. dev. 


UE noise figure 


7 dB 


System bandwidth 


10 MHz 


SNR to rate 


3 dB from Shannon, max 4.8 bps/Hz |35| 


Antenna Pattern 


macro: 3 GPP standard model [23] 
pico: omnidirectional 


Traffic model 


Full buffer 



C. User Association Optimization 

As discussed above, the single path constraint $FI\ is non- 
linear and cannot be placed into the framework of the general 
flow and interference control problem in Section [TTJ However, 
a standard method |T2|-fT4| is to solve the multiflow upper 
bound and then simply the "truncate" the solution to remove 
all but the best path. In most cases, the multiflow optimization 
yields solutions that are zero in all but one path (the book |T4| 
provides conditions under which this property is provably true) 
and therefore the truncation often does not significantly alter 
the optimization result. A similar truncation method is used in 

V. Numerical Simulations 
A. Lightly Loading Heterogeneous Network 

To demonstrate the methodology we simulated the algorithm 
in a two-tier heterogeneous network with macro- and picocells 
loosely following the 3GPP evaluation methodology in (23), 
specifically the macro and outdoor remote radio head/hotzone 
scenario in Configuration #1. The parameters are shown in 
Table [I] The main difference is that we consider a much greater 
density of picocells (10 instead of 4 picos per macrocell), 
making the number of UEs per cell small. We will see that 
under this lightly-loaded scenario, the power optimization can 
offer significant benefits. 

We focus on the problem of optimized user association (as 
opposed to multiflow) where each UE must select exactly 
one cell - either macro or pico - since this is main concern 
in the cellular standards. To apply the proposed optimization 
method to the user association problem, after the drop, each 
UE selected three candidate BSs: the strongest macrocell 
and the two strongest picocells. The algorithm was then 
run to maximize a proportional fair utility, Ui(r) = log(r), 
which corresponds to the sum-log rate. Fig. [3] compares the 
resulting rate distribution from the optimization algorithm 
against simpler, but standard, cell selection policies. The curve 
labeled "optimizer" is the optimization run, but keeping all 
the cells transmitting at maximum power. In this case, the 
spectral efficiency is constant and the algorithm reduces to 




10"' 10" 

Throughput [Mbps] 



Fig. 3: Rate distribution comparing various user association 
algorithms for 10 UEs per macrocell. Other parameters given 
in Table |I| 



the method of |4]. The curve labeled "optimizer+PR" adds 
power reduction where reducing the bandwidth allocated in 
a cell causes its power to reduce. For comparison, the figure 
also plots the "macro only" and "pico only" where the UEs 
use only the macro- or picocells. Both of these options give 
very poor rates. The range expansion methods (labeled RE) 
select the strongest cell - pico or macro - but the picocells are 
given a fixed bias in the received signal strength to encourage 
mobiles to select the smaller cell. Manual trials with different 
bias levels found an an optimal bias of approximately 6 dB. 
Fig. [3] shows the rate distribution of RE with bias of and 6 
dB. 

Table [II] shows the numerical values from the simulation. 
With the light loading of 10 UEs per macrocell (the case 
plotted in Fig. [3] we see that the proposed optimization with 
power reduction offers an approximately 23% increase in cell 
throughput and 42% increase in the 5% cell edge capacity. 
While not dramatic, the gains are larger than previous studies 
(3), (4) which showed little gain over simple RE with the 
correct bias. The reason for the larger gain here is that, since 
we are considering a low UE to cell ratio, some picocells can 
be chosen to be unused and can reduce their power. Indeed, if 
we go to a higher loading of 25 UEs per macrocell (2.5 UEs 
per picocell), the gains from power reduction are reduced since 
most of the picocell become fully utilized. 

B. Backhaul Network Constraints 

In addition to the interference control, a second appealing 
feature of the proposed user association methodology is that it 
can account for backhaul network constraints. Such backhaul 
constraints may arise when open femtocells are used for cov- 
erage. Femtocells are not connected directly into the operator 
core network, but instead use third party backhaul which may 
be variable in quality. To simulate this scenario, we employed 
the same assumptions as before, but we placed a 1 Mbps limit 
on 50% of the picocell links (the picocells now modeling large 
numbers of open femtocells). Fig. [4] shows that the gains from 
optimized user association in this scenario are even larger. 



TABLE II: Simulation Results 



Method 



10UE 



25UE 



Sp. Eff. 1 



Rate 2 



Sp. Eff. 



RE [6dB] 



0.608 



0.608 



0.909 



Opt. 



0.625 (2.8%) J 



0.745 (22.6%) 



0.940 (3.4%) 



Opt. + PR 



0.748 (23%) 4 



0.866 (42.4%) 



0.989 (8.7%) 



5% Rate | 0.351 | 0.418 (19.1%) 
1 : Spectral efficiency in bit/sec/Hz 
2:5% Cell-edge rate in Mbps 
3,4 : gains compared to range expansion [6 dB] 



Max-RSSI vs. Range Expansion vs. Optimizer (P=0.5) 



0.435 (23.9%) 




10- io u 

Throughput [Mbps] 



Fig. 4: Rate distribution of various user association algorithms 
with backhaul rate limits on 50% of the picocells. 



Indeed, for 10 UEs per cell, there is an 28% increase in average 
user throughput and 112% increase in the 5% cell-edge rate. 
In addition, Fig. [4] shows that RE gives a minimal gains in 
this heterogeneous network condition since RE is a power- 
based selection scheme and does not account for the backhaul 
constraints. 

VI. Conclusion 

We have presented a general methodology for joint inter- 
ference coordination and user association problems in cellular 
networks. The methodology is extremely general in that it 
can incorporate interference constraints described by an lin- 
ear mixing model followed by an arbitrary interference-to- 
rate mapping. Although we have only considered bandwidth 
allocations in this paper, based on the results in |T5| , we 
believe that we can extend these results to more complex IOC 
schemes as well such as subband allocations and beamforming. 
Also, while we cannot find a provably optimal solution, the 
dual decomposition method enables efficient implementation 
of suboptimal solutions via augmented Lagrangian techniques 
as well as computation of upper bounds on the maximum 
utility. 
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